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ABSTRACT 


It is proposed to explore the feasibility of using 
24 -hour, elliptic inclined orbits for satellite oonmunications. 

The study has been initiated with an analysis of the effect of 
orbital geometry on the various important aspects related to comiau- 
nications satellites, e.g., boosting energy req.uirements, visibility 
and changes in peak power and total energy required for comiminioa- 
tions. The analysis shows that the use of the proposed non-stationery 
orbit can lead to significant fuel savings. This may hov/ever result 
in some increase in communication power requirements. 

The major challenge in the use of the proposed concept 
appears to be in designing a suitable passive attitude control scheme 
which can provide for the changing satellite attitude so as to 
achieve its fixed orientation with respect to the moving line of 
sight. Fortunately, it has been possible to show that the solar 
radiation prossuro can be effectively employed for achieving this 
objective to a high degree of pointing accuracy. 
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1 . nWROIiUCTION 


'Vhat we are building now is the nervous system of mankind, 
which will link together the whole human race for better or worse, in 
a unity which no earlier age could have imagined... " said Arthur C. 
Clarke who first conceived of communication satellites more than 3^ 
years ago. 

i 

What is a communication satellite ? The communication satellite 
can be looked upon as a space relay moving in orbit, receiving infor- 
mation from a ground station and transmitting it to another ground 
station. In several respects, space communication is simpler than 
conventional communication which is severely handicapped due to the 
problems like tropospheric scatter, ionospheric radio or under water 
acoustic communications. As a matter of fact, today - 30 years after 
the prophecy by Clarke - the satellites appear to have already become 
an important link in global communication net'work. 

The various satellite systems which have been proposed for 
communications are % 

(i) Passive reflector satellites, 

(ii) Active satellites in low and medium altitude orbits, and 

(iii) Active satellites in a 24 -hour synchronous orbits. 

The merits and demerits of the various satellite communication 

2-4 

system have been examined by several investigators . Of the Several 
concepts proposed for space communication, the one most commonly used 
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and now almost nniversally accepted requires placing the satellites 
into perfect geostationary synchronous orhits. The single over-riding 
consideration that seems to have weighed heavily in its favour is 
that satellites that appear absolutely stationary in the sky from the 
earth allow the use of simpler, less expensive ground stations in 
communications network. However, there are some other important 
considerations which appear to have been ignored- This may perhaps 
be attributed to the undue emphasis being laid on the aspect of 
operational simplicity. The technological problems associated with 
the non-st at ionary satellite alternatives and the likely difficulties 
of excessive demand on altitude control in ecliptic orbits may have 
been other possible factors. 

The synchronous geostationary satellites although ideally 
suited to uninterrupted communications are difficult to achieve and 
maintain particularly over the extended periods of operations. This 
difficulty arises due to several reasons : 

(i) Large orbital deviations are caused due to the various 

perturbing forces. The frequent orbital corrections are 
therefore needed leading to rather stringent station-keeping 
requirements if the full potentialities of the geostationary 
satellites are to be realized. 

(ii) To inject the satellites into the geostationary orbit, 

large multistage rockets with high boosting capability are 
needed. The use of a non-equatorial launching site 
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(unavoida’ble if the launching has to he done from within India) leads 
to further increase in expenditure of energy, the additional energy 
being required for affecting the necessary final transfer maneuver for 
orienting the satellite orbit to the equatorial plane. This further 
affects the overall economics of the geostationary satellite adversely. 

5 

It may bo pointed out that a recent study by Dial and Cooley 
has established that the net on-station weight capability of satellite can 
be increased considerably by using a 24-hour, inclined, elliptic orbit. 
This analysis is of considerable significance as it suggests the 
possibility of affecting major fuel saving through the use of inclined 
elliptic but 24 -hour orbit for satellite communication. 

Hero, it is proposed to explore the alternative possibility of 
using elliptic, non-equatorial, 24-hour satellites for uninterrupted 
communications. It is felt that this may in general lead to conside- 
rable saving in terms of booster power requirements without compromi- 
sing to any appreciable extent the communication performance chara- 
cteristics. 

It is felt that the apparently simpler and economical proposition 

if technically feasible would open up a multitude of possibilities 

which are particularly significant in the light of growing problem 

6 

related to the overcrowding of the limited equatorial space with 

7 

satellites and of course the associated poiitical/legal implications. 

Hence, it is proposed to study some important aspects related to the 
use of 24 -hour inclined, elliptic orbits for satellite communications. 
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Chapter 2 Gxaminos effect of orbital geometry on the requirements of 
booster capability, communication power and energy, and -visibility of 
the satellites from ground station. 

It may, ho’jever, be pointed out that one of the most challcng- 
ing problems in the use of the proposed orbits for communications is 
likely to be the control of attitude motion due to excitation induced 
by eccentricity. This problem is further complicated by continual 
periodic drift of the satellites towards or a.way from local -vertical 
through the ground station. This brings us to the important problem 
of control of satellite orientation along the line of sight i.e, the 
line joining the satellite and the ground station executing periodic 
oscillations about the local vertical. This phase of investigation 
aims at exploring the possibility of using solar radiation pressure 
for achieving this variable attitude control. The problem has been 
studied in order of increasing complexity. 

In Chapter 5> o,n attempt has been made to synthesize a suitable 
but simple control law -iirhich can enable the proposed solar controller 
to orient the satellite along the local vertical with a high degree 
of pointing accuracy. 

Chapter 4 deals with the more complicated problem of developing 
a suitable, simple solar controller model for achieving the variable 
attitude control of satellite along the moving line of sight. 



2. EFFECT OF OEBITAL GilOMSTRY OH EKBR&Y AND 
GOI&IimiCATIQHS REQ.UIBEMKMTR . 


The various aspects of satellite motion likely to affect the 
accoptahility of its orhit for communication are 

(i) boost- energy requirements, 

(ii) visibility of satellite from ground station, and 

(iii) communication power and energy requirements. 

The study has therefore been initiated with an analysis of the 
effect of orbital gooim^try on those factors. 

2.1. Boost -Energy Roquiremonts ; 

The analysis compares the effect of orbital eccentricity and 
inclination on boost-energy requirements for the various 24-hour 
orbits. Since injection of satellite into the final desired orbit 
is in general preceded by establishing the circular parking orbit 
first, attention has been focused on computing the total velocity 
increment required in transfer of satellite from the parking orbit 
to final desired orbit. Since the objective of analysis is limited 
to only a prcliminaiy estimate of impulse requirement and their 
comparison, the 2 -impulse, minimum energy i.e. Hohmann transfer 
trajectory has been assumed for transfer to the perfect geosynchronous 
orbit. Furthermore, the final impulse at the apogee of the transfer 
orbit is modified to provide for the change in inclination of the 
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orbital plane. In the alternative case of transfer to the general 
24-hour elliptic orbit, a single impulse is assumed to act at the 
perigee (Pig. 2.1) of the desired elliptic orbit. 

The dimensionless expressions for the total velocity 

increment required in transfer from the parking orbit to the final 

circular or elliptic orbits with common 24-hour period are found 
* 

to be j 


AYq = /y /a^ C 2 n 2 - 2112/(1+ n 2 ) - ^/n^ 



Using these expressions, it is possible to compute the percentage 
saving involved in the use of the proposed elliptic orbits instead of 
the perfect circular orbit. These results are shown in Table 1. 

It is thus apparent from these results that the use of elliptic 
inclined orbits can enable a significant saving in booster power 
requirements. The saving increases v/ith increase in eccentricity as 
well as inclination. 


* 


See the Appendix 
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2.2. Visibility of Satellites from Ground Station . 

An important condition for use of satellites as communication 
relay without time delay ia its yisibility from the ground stations 
being used for the purpose. This imposes constraints on the maximum 
permissible orbital eccentricities that can be used with satellites 
in 24 -hour orbits. An attempt has, therefore, been made to obtain a 
theoretical estimate of this upper limit constraint on eccentricity. 
By analyzing the geometry of motion (Pig. 2.2) it is easy to see 
that the yisibility of satellite from the ground station is assured 
provided the folloviring inequality constraint is satisfied 

cos 0 > r^/r ... ( 2 . 3 ) 

where, cos 0 = r . r^ /( |r | jr^l ) ... ( 2 . 4 ) 

On evaluating the expressions for cos 0 in terms of the 

8 

orbital elements, this inequality takes the form 

r pt COS u cos(n-6 ) - sin u cos i sin ( fi - d)} cos <j) 

+ ( sin u sin i) sin ij* ]] > 1 ...( 2 . 5 ) 

where r^^ = 

XX = ( 1 ) 4* 0 

6 = 60 + tOe t 

«o = 


longitude denoting the initial position of the ground 
station with respect to vernal equinox 
latitude of the ground station 



“C- e sin 6 /l-e^/ (i+e cose) + 2 tan’’V (l-e)/(l+e) tan 8 /2} J 

••• (2»4) 

This ooustrEitii; ceu also be written as 

+0— 5 q) cos <(» + { (l- cos i) sin (fi- 6^ - a)gt)sinn} 
cosi^+ sin n sin i sin > 1 ... (2.7) 

It is apparent that the visibility of the satellite would be 
influenced by it,s initial longitudinal position Tirith respect to that 

of the gioxmd station. Normally, the relative position of satellite 
injection into the final orbit would be so chosen as to ensure the best 
visibility conditions, leading to a considerable simplification ill the 
above relation 

^n [pos(6 - “e't) cos - (l_ cos i) cos ()) - | sin i sin <(. ( J > 1 

- • • • (2 *8^ 

Using the expression for obtained earlier, it can be shown that 


cos 


( 6 - ajgt) 1 - 4-e^ sin^e > - Ae^ 


... (2 .9) 


Substituting (2.9) into (2.8), the inequality takes the 


form 


/ 


1 - 4e > 1 - cos i + ( sin i tan 4) | + sectf) / (a^(l-e)} 


...( 2 . 1 ») 


where 




*= ^ *628 
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The result is of considerahle general significance hoing 
independent of the several orbital elements, e.g. 0 , oi > > 6o * 

Furthermore, it provides useful design information giving the maximum 
permissible eccentricities, consistent \7ith visibility, as a function 
of orbital inclination v/ith equatorial plane as well as the latitude 
of the ground station. Results are presented in Table 2. 

From these results, it is apparent that the use of largo 
eccentricities is possible for communications even for relatively 
large inclinations and for large ground station latitudes. For 
example, it is possible to ensure visibility for eccentricities upto 
0.55 oven in cases whore i and <i> become as large as 30°. 


2.5. Gommunication PoTrer Requirements s 

Recognizing that the communication povjer required for trans- 
mission from or to satellite is proportional to the square of the 
distance between tho satellite and the ground station, it can be 
shown that 


P 


max 





... (2.11) 


where 


P == Communication Power required for geosynchronous orbit, 
ref 

P = TTrnm Commiinication Power required in case of alternative 

max 

elliptic, 24-liour orbit. 

This expression can be rewritten in a more convenient form as > 



MAXH-IDM FBEUISSIEffi ECGEKTRIGITIbS FOR Min'i'EREUPi'ED VISIBILITY 
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W^ref - D + 0/(1- ... (2.12) 

or 

W^xof - [3 + 1 .176 K ^ ... ( 2 . 12 a) 

The results obtained using this relation are sho?m in 
Table 5» It may be noted that the comiminica.tion power increases 
with increase in eccentricity. However, this increase in peak 
pOYiTor requirement may not be able to offset the gains of fuel 
saving as it would also be accompanied by the corresponding 
decrease in the pov/er required at the points of closest approach. 
A more realistic comparison would, therefore, bo provided by 
estimate of the total communications energy required* 

2 .4 ■ Communication Energy Requirements 

A comparison of energy requirements for communications for 
the various eccentricities can be made using the relations 

21T 

E a / P de ... (2.15) 

0 

where the power P can bo represented by the relation 

P a (r - ••• (2.14) 

Substituting (2.1 4) into (2 .15) and carrying out the necessary 
manipulation, we get 
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Ea J (r^_ 2r rg + r|)d0 (2.I5) 

0 

On evaluating the integrals appearing in this expression^ it is easy 
to show that for elliptic orhits 

Ee a 2Tr [a^ _ 2a rg A- e2 + D (2.16) 

Por circular orhitsji i.e, when e = 0, this expression reduces to, 

Ec« [a^ _ 2a rg + ^ (2.17) 

On. dividing (2.16) Toy (2.I7) vre get 

Eg/Eg = + ( 1 - A- e 2 )/(a^-l )2 (2.13) 

or 

Eg/Eg = 4 -e^ + 0.032(1- M-e^) 

Tafele 4 shov/s the effect of orbital eccentricity on increase 
in oonunanication energy in elliptic orbits. It may bo noted that 
the total energy requirement shows a slight decrease with increase 
in eccentricity. The gains offered by this marginal decrease in 
communication energy are however not very significant. The requirement 
of increased peak power for communication does appear to an additional 
problem. However if the uninterrupted communications is not essential 
there is a possibility of combining the satellite communication 
during the periods v/hen the satellite is relatively closer to the earth 
with other mission objectitres. In such a situatiouf the main 
advantages of elliptic orbits namely the higher peak power require- 
ments for oomiminication can be substantially overcome. 



ATTITUDE COMTROL OP SATELLITES IN ELLIPTIC ORBIT S USING 
SOLiR RADIATIoFmSSURE 


5.1 Introduction 

For many space applications such, as communications, weather- 
forecasting, survey of earth-resources, it is necessary to maintain 
the satellites at a fixed orientation with respect to the earth. 

The use of active attitude control methods although desirable for 
achieving high degree of pointing accuracy leads to penalty in terms 
of increased vjeight, reduced space- availability on board and limited 
satellite life-time . That perhaps explains why there have been 
numerous attempts at developing attitude control methods utilizing 
environmental forces . 

The possibility of using Solar Radiation Pressure (SRP) for 

9-14 

attitude control has been explored by several investigators . Even 
though effectiveness of this concept has been clearly established, 
particularly for satellites in high altitude orbits, it appears to 

have found only limited acceptance by the spacecraft designers. This 
lack of enthusiasm for the solar attitude controllers may be partly 
attributed to the availability of the well proven conventional passive 
and semi-passive alternatives using gravity gradient, satellite spin, 
magnetic tcrduiag etc. The oompleri-ty of the proposed controller models 
and their operation has perhaps been a major detrimental factor as 
well. It seems that from the very beginning, there has been nndne 



16 


emphasis on using the combination of proportional and derivative 
feed-back control policies for achieving high degree pointing 
aQ^curacy . In the process, certain obvious possibilities emphasizing 
the simplicity of control and operation appear to have been overlooked. 

This chapter investigates the feasibility of developing a 
simple solar attitude controller for satellites operating in elliptic 
orbits. Attention is focused on evolving a simple control strategy 
which can counter the adverse effect of eccentricity normally 
responsible for worsening of the attitude characteristics of the 
conventional methods. 

5*2 Formulation of the Problem : 

The formulation begins with the analysis of motion of 
a cylindrical satellite in elliptic orbit, in the ecliptic plane 
about the oarth-cGntcr 0 (Fig. ^be satellite model considered 

here carries two sets of plates made of light, rigid but highly 
reflective material (o.g. aluminized mylar membrane). The two sets 
are assumed to be mounted with their surfaces normal to the orbital 
plane. Denoted by x,y,z are the principal coordinate axes of the 
satellite with origin at the mass-centre s. The angle made by the 
z~axis with the local vertical defines librational angle, Tp . Using 
Lagrangian formulation, the governing equations of the planar libra- 
tions can bo written as 

(1 + o cos e) Ip" - 2e(l + ,p’)s±n 0+ 3 K. sin 4, cos ip =(R^/vl)q 

... (5.1) 
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where Q,^ represents the generalized foixse due to the SEP. 

For specular reflection assumed and using the principle of virtual 

work, the expression obtained for the generalized force takes the 
form 

% = (s/c' ) Ql +P^) «>^A^sin( 0 + i|^ _ 4)) I sin (0 _<j) ) I 

- (1 + P 2 ) ^2'^2 ^ +'P - <)>) loos (0+i|)_<f))| 3 

... (5.2) 

substituting the above expression for in (5.1) transforms the 
equation of motion to 

(1 + e cos 0) 4)" - 2o (1 + 4»') sin0+ 3 sin4)oos4' = 

3 ^ 

Ul + e)7(l + e cos 0)^}[C^ sin (d + ip -(p ) [sin (0 -.<!> )[ 

- C2 cos ( 0 + ill- <}> ) 1 cos ( 0 + il' - 4 > ) 1 U ... (3,3) 

whore , Cg referred to as the solar parameters in the analysis, 
are given by 

o A. 1.(1 + pp r5 j = 1,2 ...(3.4) 

3.3* Effect of Solar Radiation Pressure on Attitude Dynamics . 

The governing nonlinear, non -autonomous differential 
Equation ( 5 . 5 ) does not have a closed-form solution, numerical 
techniques were therefore used to evaluate the librational 
performance of the satellites in various practical situations. 
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Numerical integration. V'/as performed using Adams Extrapolation Method^ ^ 
in conjunction with a suitable step-size of 2". 

Figure 3*2 shows the librational response in tuo typical 
situations ^/ith different orbital eccentricities. To bring out the 
influence of the SEP clearly, the response for = C 2 = 0, i.e., 
when the solar effect is ignored, has also been included for 
comparison. It may bo observed that the SEP increases the amplitude 
of librations when is positive. On the other hand when negative 
values arc chosen for this parameter, the SEP seems to influence the 
system favourably reducing the maximum amplitude of librational 
motion. 

To understand the dynamical behaviour of satellites over a 
v;idu range of system paramotors, numerous response plots were 
obtained. The resulting information is condensed in the form of 
system plots (Pig. 3*5 - 3*5)* 

Figure 3.3 sho\;s the effect of eccentricity on the librati- 
onal amplitude in throe different situations. It may be observed 
that in general the amplitude continually increases with increase in 
eccentricity. Hero, solar parameters , Cg and the solar position 
angle <t> appear to bo the other significant parameters of the system 
affecting the librational motion. 

Figure 34 indicates the effect of the solar parameters 
, C 2 on the librational response. It may be pointed out that in 
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casG of circular orbits and for =02 = 0, the satellites with 
thoir long axes initially oriented along the local vertical do not 
execute the librational motion unless disturbed. In contrast the 
response results for any other non— zero combination of the solar 
parameters exhibit considerable amplitude of librations thus estab- 
lishing the detrimental influence of the SRP in the circular orbits. 
However, it is interesting to note that for elliptic orbits, the 
minimum amplitude does not correspond to the case = C 2 = 0, i.e., 
when the solar efCoct is ignored. Furthermore, through a judicious 
choice of the solar parameters, it is possible to reduce the libra- 
tional amplitude substantially. The significance of this result 
cannot bo over-emphasized as it suggests an attractive possibility 
of countering the adverse influence of eccentricity on attitude motion 
of satellites. 

The effect of the varying solar position angle is indicated 
in Figure 5*5* noted that the solar position has conside- 

rable influence on the librational amplitude of satellites, 
particularly in elliptic orbits. This suggests that a change in <p 
signifio.-mtly affects the choice of the solar parameters for achieving 
an improvement in the satellite librational performance. 

Thus, the numerical results presented here clearly establish 
the possibility of a.chieving considerable reduction of the maximum 
librational amplitude through on appropriate choice of the solar 
parameters. This beneficial aspect of the SRP appears to have been 
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cither missed or entirely ignored so far. On the other hand, an 
impression scums to have gained ground that the uncontrolled solar 
effect is always detrimental to the attitude performance. This may 
perhaps partly explain why, in the recent developments of the various 
solar controller models, the emphasis has "been essentially on 
regulating the solar torque through the moving surfaces using complex 
feedback control systems. 

The system plots obtained here can no doubt be used for 
selecting the appropriate values of the solar parameters in several 
cases. However, to numerically generate the extensive data to enable 
the choice of those parameters covering all possible situations is 
neither fcjasiblo nor desirable. So, it was now proposed to treat the 
problem analytically T/ith a view to develop the suitable criteria 
facilitating a judicious choice of the solar parameters. 

5.4 Synthesis of Control Strategy using Auuroximatc Inalytical Approach . 

It was felt that in view of the complexity introduced by the 
nonlinear and non-autonomous nature of the problem, it would be 
essential to make some reasonable assumptions for treating the 
equations of motion analytically. It is easy to see that for the 
practical situation of small amplitude librations assumed here, the 

16 

governing equation can be written as 

(1+ e cos e)f - 2gV sine+ 3 K^ sin cos 1(1= 2 e sine + 

{ (1+0) ^/(l+e cose)^} [c^ sin( |sin(e - <J) ) | 

- G 2 cos( e ) lcos( e - <}> ) 1 H 


... (5.5) 
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It may ‘bo pointod out that in. general the deterioration of 
the attitude control performance characteristics is mainly due to the 
periodic excitation represented hy 2e sin 9 . Attempt was therefore 
made to explore the possibility of getting rid of this term with the 
help of the l.ast two terms representing the solar pressure effects. 

To achieve this objective, the restnmoturing of the solar teims was 
considered so that \?ith a suitable choice of C..J, the right hand 
side in the equation could be forced to disappear. Fortunately, 
Fourier expansion of the terms appearing v;ithin the square brackets 
in Equation (3.5) suggests that they can be effectively described by 

16 

their first harmonic . On incorporating this approximation, the 
Equation (5»3) simplifies to 

(1+ o cos 6)(|) - 2ei()' sin9 + sini|)cos4’= 2e sin 6+ (s/jir) 

{(l+e)/(l+o cose) [~C..| sin (9-iJ>) |sin (0 -({>)] 

_ C2 cos ( e _<{) ) |cos (e - <}) ) I ]] ... (3.6) 

On separately grouping the terms involving sin 0 and cose , it was 
found thiit the right hand side in the Equation (3*5) vanishes when 

= -0.75'^ {g/(1+ e)^'} (1+e cose)^ cos (|) 

C2 = +0.75Tr {o/(l+e)^ } (l+c cos 9)^ sin ()> ... (3*7) 

These results are of considerable significance, as they 
clearly bring out the effects of the orbital eccentricity as well as 
the sun— position angle on the desired instantaneous values of the 
solar parameters. It may be pointed out that by suitably regulating 
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the movement of the pairs of the solar surfaces and 

respectively along the z and x-axes (^ig. 5*1)> > ^2 hence 

, O 2 can he varied as desired. Thus using these relations, it 
may he possible to control the solar parameters suitably in order 
to overcome the adverse effect of eccentricity. 

A careful examination of the control relations reveals that 
the changing satellite position in orbit is likely to be of little 
consequence to the selection of , Cg, particularly for low values 
of orbital eccentricities. Hence, these relations can be considerably 
simplified 5 

= -0,75ir {e/(H-e)^ } cos <> 

Cg = + 0 . 75 TT {e/(l+e)^ } sin (J) ...( 5 . 8 ) 

To examine the validity of the approximate analytical 
approach, the above ejcprossions developed for were substituted 

in ( 5 •5) and the resulting equation of librational motion was 
integrated numerically. Some typical response results thus obtained 
are shown in Figure 5*6. It may be observed that the control schemes 
developed analytically are quite effective in limiting the amplitude 
of librations to much Icwjor values. Even v^ith the highly simplified 
control of solar parameters represented by (5o8) where = C^(<j>), 

C 2 = CeC^) j the amplitude reduction is found to bo considerable. 

The more stringent control scheme represented by (5*7)» whore 
C-i =C^(0,<j)), C 2 = C 2 (ej({)) further improves the controller 


effectiveness . 
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To assess the influence of the important system parameters 
on controller effectiveness, the response results 'vfere analysed in 
details. Piguro 5*7 shows the effect of the solar position angle on 
the system "behaviour. As the position of the sun changes continuously 
in the ecliptic plane, this study would bo particularly useful in 
predicting the long range satellite librational dyncmics. The analysis 
shows that the changing position of the sun does not significantly 
alter the attitude control characteristics, pa.rticularly for low 
values of occontricities. 

Figure 5 -8 exhibits the effect of orbital eccentricity on 
maximum librational amplitude. It may be observed that both the 
control strategies are quite effective in reducing the amplitude of 
motion. For continuously increasing values of eccentricities, the 
effectiveness of the control schemes progressively declines. The 
degradation in performance is particularly severe Ti/ith the simpler 
control scheme. This limits the applicability of the 4 * -sensitive 
control to only slightly eccentric orbits. However, usefulness of 
the proposed, control scheme sensitive to both 0 as well as (}i , 
extends to much larger range of eccentricities. It is interesting to 
note that even for eccentricity as largo as O. 5 , it can limit the ! 

amplitude to well Virithin 5"« viev; of the recently proposed use of 

highly elliptic geosynchronous orbits in order to maximize the net 

5 

on-station v;oight capability , this study may be particularly significant 
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Attention was now focused on the comparative assessment of 
the overall attitude performance characteristics associated with the 
two proposed control laws for the particular case of satellites in 
near-circular orhits. It was felt that, in this situation, even 
though the control operation sensitive to hoth 0 and provides 
slightly "better pointing accuracy, the tj) -based action ¥/ould result 
in considerably lower demand on the control system. It may be pointed 
out that the simple, <i> -sensitive control involving the single 
variable is easier to realize than the other one dependent on both 
6 and <l> . Furthermore, as the period of terns such as cos and 
sin (j> is one year -much higher than the orbital period of the 0-terms 
involved in the control relations, it is apparent that another major 
advantage of the simple 4 -sensitive control function lies in its 
much lower frequency. It may be possible to use this fact to 
advantage by effectively replacing the continuous periodic control 
action by a limited number of presfecified discrete operations over 
each cycle. To assess the effectiveness of this approach, the 
response plots were obtained for the varying degrees of discretization 
of the proposed ^-sensitive control criteria. The results of this 
analysis are shown in Figure 3 •9* Here, represents the solar 

position angle with respect to vdiich the solar control surfaces are 

t 

assumed to be set initially and A4> denotes the maximum permissible 
deviation in before permitting the resetting of the solar surfaces. 
It may bo observed that the discretization leads to considerable 
deterioration in the attitude control characteristics. However, for 
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values of &<}> representing small changes in ‘f’ of say upto + 5 '’» "the 
minor degradation in librational performance is likely to be •well 
within the tolerance limits for many space applications. Ihis result 
is of considerable practical significance, as it establishes the 
feasibility of employing the solar surfaces suitably mounted on 
satellites to achieve substantial reduction in librational amplitude, 
even -without exercising control during the long intervals extending 
upto as much as a T/eek. The proposed control approach thus suggests a 
simple method for testing the concept of achieving libra-tional 
control using the SEP with minimal changes in satellite design. 

5.5 Effect of Generalizing the Control Law on Satellite Attitude 
Performance . 

The open-loop control policies suggested here altho-ugh quite 
effective in limiting the librational amplitude to low values, does 
not provide for damping the excited librations. It is, therefore, 
proposed to combine them -with the proportional and derivative control. 
The resulting control laT<\r can be -written as follovj-ss 

•where 

or C^( 6 ,((>) 

*^21 ~ 02(9, <l>) 

I 2 kiT< 0 - (J ^ ( 2 k+l) tt 

+ (vi^ V+ I ( 2 k + 1 ) -rr < 0 + ip -•> <t> 1. (2k+2) 
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+ (^2’^'+ '^ 2 '^) 5 (2k- i)iT< 0+t}) _(J l(2kH-|-) TT 
- (^ 2 '^'+ ^ 2 '^)°, (2k+-|)Tr < 6+i{;_(|> ^ (2k+-|> 


... (5*10) 

with 1 0 ^ 2 ! — ^ "being an integer. 

The various results are obtained hy solving Equation (5*9) » 
using the earlier numerical techniques in conduction vrith a suitable 
step size of 0.5° . 

Figure 5. 10 shoYirs the librational response of satellite in 
some typical situations. It may be noted that using the control relation 
Cj = Cj( <{> ,^,’1'' ) obtained by the combination of (5.^^) and (5 .10), loads 
to significant reduction in the steady state librational amplitude; 
while (5 *9) reduces the librational amplitude even further. 


To understand the behaviour of system and control parameters 
on controller effectiveness, numerous response plots were obtained. 
The resulting information is condensed in the form of system plots 
(Pig. 5.11 - 5.14). 


Figure ^,11 shows the effect of mass distribution on the steady 

state librational amplitude in several situations. It may be observed 

that the inertia parameter, has considerable influence on the 

amplitude for = 0, i.e., when only open-loop control 

C - = C.(‘i>) or C. = C.(6,<{>) is employed. For the feedback control 
0 J J J 

represented by nonzero values of y. and v., however, the chaaoge in K. 

It] 3- 
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appears to have no offoct on the systora plots. Furthermore introduc- 
ing the additional terms to offset the effect of excitation induced 
hy eccentricity reduces the lihrational anplitude considerably. 

In general, the amplitude of oscillations is found to be lowest when 
the more stringent control strategy repr^-sented by C. = C . (® ) 

is employed. 

The influence of orbital eccentricity on the amplitude of 
motion is shown in Fig. ^,12. It may be noted that in general, 
increase in the eccentricity adversely affects the lihrational 
performance . However, the suggested modifications in the control 
strategy substantially improve the controller effectiveness and limit 
the amplitude to much lower values even in case of relatively large 
eccentrities. 


Figure 3*13 shows the effect of control gains vi. and v. on 
lihrational amplitude. It may be observed that a judicious choice 
of the gains can minimize the attitude variations considerably. Even 
for the best choice of controller gains, the suggested control 
modifications appear to bring about a considerable reduction in the 
amplitude of lihrational notion. It is apparent that even for 
eccentricities as large as 0.1, the proposed controller can limit the 


amplitude approximately to within 0.1°. 


Figure 3. 14 shows the effect of C. and the initial 
impulsive disturbance general an increase in C . affects 

the attitude performance favourably to a point beyond T/hich the 
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dependoncG 1200011103 rathor weak. It may to otsorved that 7/ith the 

introduction of the suggested modification in the control law, it is 

possihlo to use much smaller -values of C . 

;]max 

The proposed controller presents an interesting possibility 
of controlling satellite orientation in elliptic orbits. The reali- 
zation of the control is easy and the size of the plates required quite 
modest. A preliminary estimate using the data of HTTELSAT IT series 
of satellites shows that the plate size of O.4 vfith permissible move- 
ment of 20 cm. is adequate for orbital eccentricities of 0 . 1 . For 
larger eccentricities, this requirement increases in nearly the same 
proportion. The semi-passive character of the system promises increased 
life-time with overall reduction in the cost. 

3.6. Concluding Remarks . 

The important features of the analysis and the conclusion 
based on them mo-y be summarized as follows 

( i ) Contrary to as is generally believed, the SRP need not be 

al-ways detrimental to the satellite attitude response. In fact, 
there are several situations, iThore for a proper choice of the 
solar parameters, the radiation pressure can substantially 
reduce the maximum librational amplitude of satellites in 
elliptic orbits. 

(ii) The system plots, developed numerically, summarize the satellite 
librational response characteristics over a -ivide range of 
parameters. Through these plots, it is possible to select the 
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appropriate values of the solar parameters minimizing the 
liTorational amplitude in several oases. 

(iii) The approximate analytical approach adopted here seeks to make 
use of the solar terms in the equation of lihrational motion 
to eliminate the major excitation effect of eccentricity 
normally responsible for the worsening of the attitude control 
characteristics of the conventional methods. The two control 
strategies thus evolved are found to be quite effective in 
restraining the satellite oscillations. 

(iv) The proposed control schemes are much simpler as they are 
dependent only on 0 and ip , which in turn can be expressed as 
a function of time alone. This makes them p sarticularly suited 
for preprogramming of the control function unlike the earlier 
complex feedback control approaches sensitive to the libra.tional 
disturbances and requiring on-board sensing and computation 

of the error signal for actuating the controls. 

(v) Using the rele-tively more stringent control law governed by 
the instantaneous values of the position of the sun (<J)) and 
satellite (6) leads to drastic reduction in maximum librational 
amplitude. It considerably extends the effective range of 
possible eccentricities for stable satellite operation. It is 
interesting to find that even for eccentricity as large as 0.5, 
the proposed controller can limit the maximum amplitude to 
well within 5°. This suggests the possibility of using it as 
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an auxiliary control devico for improving the overall lihrational 
performance characteristics of satellites in highly elliptic orbits* 

The proposed control mechanism may thus have far reaching implications. 
In viCY/ of its ability to overcome the adverse effect of eccentricity 
to a large extent it may lead to greater acceptance of highly 
elliptic orbits, particularly for multi-space missions. 

(vi) The simpler <!> -sensitive control strategy is q,uite effective for 

satellites in slightly eccentric orbits. Furthermore, it is possible 
to effectively replace this continuous periodic control action by a 
limited number of prospecified discrete operations ovor en.ch cycle. 

The analysis ostchlishes the feasibility of employing the suitably 
mounted solar surfaces on satellites to improve the lihrational 
performance even Yrithout exercising control for long intoirvals. The 
proposed control approach thus suggests a simple method for testing 
the concept of lihrational control using the SEP with minimal changes 
in satellite design. 

(vii) Combining the suggested open-loop control with proportional and 

derivative feedback control appears to bring abmt a considerable 
reduction in the amplitude of lihrational motion. For example* the use 
of the modified control law enables the solar controller to limit the 

O 

lihrational amplitude to within approximately 0.1 , even for eccentri- 
cities as large as 0,1. 



4* ATTITUIE COm'HOL OF SATELLITES Bf HOK-STAIIOIARY 2A-HOIIR OB3ITS 


4.1 Introduction 

The advantages of having oommanications satellites in a 
perfect geosynchronous orbit are well understood. These orbits not 
only enable the use of simple ground stations but also facilitate 
the use of relatively simpler subsystems for station-keeping, 
attitude control and communications. It is, therefore, natural to 
conclude that any major deviations from the well established practice 
of using perfect synchronization between the satellite and the ground 
station is going to pose a major challenge on several fronts and a 
major effort is called for to meet these challenges successfully. 

In particular, the problem of using elliptic, inclined, 
24-hour orbits for communications is significantly complicated by 
continual periodic drift of satellites from their synchronous position 
directly over the ground station. The non-stationary nature of the 
satellite causing periodic oscillation of the line of sight about the 
local vortical leads to a continual periodic change in its orientation 
as viewed from the ground station. 

This problem can perhaps be easily handled through the 
active attitude control method. However, on-board energy requirements 
involved here are likely to offset the advantages of economy in 
booster capacity needed for placing the payload in the desired orbit. 
It seems, however, possible to overcome this difficulty through the 
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passive control schemes using solar radiation pressure. It is, 
therefore, proposed to explore the possibility of evolving a near- 
passive solar attitude controller which can vary the satellite 
orientation continuously so as to compensate for the steady or 
periodic satellite drift relative to the ground station. 

4*2 Formulation pf the Problem : 

Geometry of the satellites assumed to be moving in 
elliptic equatorial orbit is shown in Pig. 5 .I.. The angle, 5 made 
by the z-axis with the line of sight defines satellite orientation 
as apparent from the ground station. Recognizing the relation 
beti/ireen this apparent angle, 5 and librational angle, 41 as 

tl/ = C + a ... (4.1) 

where 

a = tan”"' { sin (0 - a)gt)/(r^- cos(e- co^i)) )■ ... (4.2) 

Using ( 4 . 1 ), the Equation (5.5) can be written as 
(1+C cos )5" - 2e(l+C') sin + 5K^ sin^cos C 

= { (l+e)^/(l+e cosS)^ }[C^ sin (6+g +a ) ] sin(0+ 5 + a - 4 ) | 

- C 2 oos( 0 H- 5 + a - (}) ) I 003(0 + 5 + a -4) )| 3“ ("I + ® cos 0 )a" 

+ 2e sin 0 a' + 5K^ sin 5 cos 5 - 3\ cos ^ sin 4 ; . . . ( 4 . 5 ) 

4.5 Effect of SRP on Attitude Dynamics ; 

To understand the dynamical behaviour of satellites over a 
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wide range of system parameters, various response plots were obtained 
numerically. The resulting information is condensed in the form of 
system plots (Fig. 4.I - 4.5). 

Figure 4*'I shows the effect of eccentricity on the 
librational amplitude in three different situations. It may he noted 
that the maximum librational amplitude increases ■vfith increasing 
value of eccentricity. The effect of solar parameters and Cg on 
the librational motion is shown in Fig. 4*2. As before, it may be 
pointed out that the minimum amplitude does not correspond to the 
case when = C 2 = 0* Furthermore, through a judicious choice of 
solar parameters, the librational amplitude can bo reduced considerablj 

Figure 4*5 shows the effect of variation of solar 
position angle, <{) . It may be noted that (f> has considora-ble 
influence on librational motion. This shows that a change in (j) 
significantly affects the choice of the solar parameters to achieve 
improvement in the satellite librational performance. 

In view of the obvious limitations of numerical approach 
leading to a judicious choice of the various parameters, attempt 
was now made to use an approximate analytical approach for predicting 
the suitable choice of control parameters. 

4.4 Analytical Apnroach s i 


Due to complexity and nonlinearity of the nature of the 
problem several simplifying assumptions have been made to treat the 
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problem analytically. Pox small amplitude motioUj it is possible to 
approximate (-4.2) by the relation 

a sin (e- Wgt) (l + e cos8)/a^ ... (4.4) 

where 

a^ = a^(l - e^) - 1 

Linearizing (4 . 5 ) in con juction with (4 .4) , equation of 
motion can be written as 


(l 4 -e cos0)C - 2e sin0 C + 3 ^^ sin ^ cos g = 2e(l- jK^/a^) sin 6 
+ {(l+e)/(l+e cosS) } ^ sin(0 - ((f)-a )) 

I sin( 9 -((|)_o))|- C2 cos(6 -(4 -a )) | cos(e -(<})-a)) 

(4.5) 

Following the steps as outlined in Chapter 5, it can be 
shown that the effect of excitation due to eccentricity and gravity 
gradient moment can be counteracted by controlling the solar parameters 
according to the relations 


= _0.75Tr 0 (l-5K./a^) cos Y (l+e 00s 0 ) 

7 7 

Cg = +0.75 ire (l-5Kj_/a^) siny ^1+e cos9 )'^/(l+8) ^ 


... (4*6) 


where y = (j) - a andot ^ (e-a^t)/a^ 


These results are of considerable significance, as they 
clearly bring out the effects of the orbital eccentricity inertia 
parameter and the sun position angle on the desired instantaneous value 
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of the solar parameters. It may he noted that hy suitably moving 
the solar plates, the solar parameters can he varied in order 

to overcome the adverse effect of eccentricity. 


A careful examination of the control relations reveals 
that the changing satellite position in orhit is likely to he of 
little, consaquence to the selection of , C2 particularly for low 
values of orbital eccentricities. Hence, these relations can he 
considerably simplified; 


= -0.75TT e(l- 5 \/a^) oos<J>/(l+e)^ 
Cg = +0.75 ■'T e(l-5K^/a.j) sin<l)/(l+e)^ 


... (4.7) 


To examine the validity of the approximate analysis, 
the above expressions for 0^, C2 were substituted in ( 4 * 5 ) "the 
resulting equation was solved numerically. The results are shown 
in Fig. 4.4. It may bo noted that the control schemes developed are 
indeed quite effective. However, the approximate control of solar 
parameters according to 0^ = C^(<f>,K^) represented by (4.?) 4 as almost 
the same effoct as that due to = 0^(6,'!’ ,Kj^) denoted by (4.6). 

To understand the dynamic behaviour of system parameters 
on the controller effectiveness, the response results wore analyzed. 
Figure 4.5 shows the effect of the solar position angle on the maximum 
librational amplitude. This study may be useful in predicting the 
long range satellite dynamical behaviour. The analysis shows that 
there is not imeh change in the librational amplitude especially for 
low eccentricities* 
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Figure 4*6 shovirs 1:11© effueb of occGti'trici'fcy on. "tlio 
ninxiiKuni libriiiiional a.mpli'tud.G . 11; may to oTssorvod "thal; tofli control 

strategies are equally effective in reducing the amplitude of notion 
for the low values of eccentricities. As the value of eccentricity 
increases, the effectiveness of control schemes decreases. Hov/evor, 

it is worthv/hile to note that the amplitude can he limited well 

0 

within 5 oven when e = 0.10. 

4 *5 • 'Effect of Generalizing the Control Laar on Satellite Attitude 
Performance 

As the open loop control policies suggested above are 
not fully effective for higher values of eccentricities, it is 
proposed to combine them with the proportional and derivative control 
law. The resulting control law can be written as % 

+ C^2 = 0.(0, <!>, i>') ... (4.8) 

v/hcro 

C.,., = C.,(<I',K^) or C.,(e,<J., K^) 

= C2(<J',Kj_) or C^iQrh \) 

C ^2 = 0.^(4', V) and Cgg = are given by (5.IO) 

Since the modified governing equation of motion does not 
have closed-form solution, numerical techniques were applied to study 
the system behaviour. Hero, a stepsize of 0.5° was found to be suitable 
for numerical integration. 
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To examine the suitahility of control given by (4.8), 
some typical response results were obtained (Fig. 4.7). I'b 't'® 
observed that the solar control given by C . = is quite 

^ t3 

effective in limiting the librational amplitude .The more stringent 

control schemes represented by C. = C 4' ) combination of 

(4.7) and (5.10) and C. = C.(6, ij) combination of (4*6) 

D t3 ^ 

and (5.10) further improves the controller effectiveness. 

To assess the effect of system parameters on the steady 
state librational amplitude, numerous response results were obtained 
which arc presented hero in a rather condensed form (Fig. 4.8). 

It nay bo obsoarvod that with the proportional and derivative control 
the steady state librational amplitude is reduced considerably. 

The amplitude can bo reduced further using the modified control laws 
C- = K. ,4','}'’) and C. = C .(6 , 4’) . Furthermore, the 

results suggest that the simpler control strategy denoted by 
C. = C .(4),K. ,4), ) is quite adequate and should hence be preferred. 

4.6. Concluding Remarks. s 

Tho important features of the analysis and the conclusion 

based on them may be summarized as follov/s; 

proper choice of solar parameters can reduce tho maximum 
lihrational amplitude substantially. 

(ii) The system plots, developed numerically, summarize the 

satellite librational response characteristics over a wide 
rango of parameters. Through these plots, it is possible to 



58 


so loot “the appropriate values of the solar parameters 
ninimiziug tho librational amplitude' iu several cases* 

(iii) Tho tvo control strategies 0. = C.(ij>,K.) and C . = C,(e, ^ ,K.) 

j 3 1 3 j i 

aro found to he quite suitable in overcoming the major 
excitation effect of eccentricity and inertia parameter. 

(iv) The simpler control strategy represented hy C - = C.('f’,K.) 

j 0 ^ 

can ho employed quite effectively for limiting the 
librational amplitude Tmthout appreciably conpromising 
v/ith pointing accuracy. Tho resulting solar controller can 
also bo used as an auxiliary control device for improving 
tho overall librational performance characteristics of 
satullitos in elliptic orbits if higher degree of pointing 
accuracy is required. 

(v) Tho addition of derivative and proportional control makes 
tho controller more effective. This generalized control 
limits tho steady state librational amplitude to well within 
5° for Gccontricities even as large as 0.20. 

(vi) Tho generalised control policy indicated by ) 

appears to bo the best as it generally leads to minimum 
limit cycle anplitude of librations. 



5. CLOSIHG EEMifflKS 


5*'' Summary of tho Conclasion.3 

As already stated the study aimed at examining the -various 
aspects of motion related to the use of 24 -ho\ir, non-stationary orhits 
for satellite comrtiunioatien. Tho important conclusions hased on the 
analysis can be summarized as follows s 

(i) Tho use of the proposed non-synchronous orbits can enable 
substantial saving in the booster capability requirements. 

(ii) Even for relatively large eccentricities, and with moderate 
orbital inclinations and ground station latitudes, the 
uninterrupted visibility is possible. 

(iii) The requirements of peak po\Ter and total energy for 
cominunications show a continuous increase with increasing 
orbital eccentricity. However, this may not pose a major 
limitation, particularly with proper scheduling of 
cominanion,tion periods and when this main objective 
satellite communication objective is combined with various 
other important missions, e.g., remote sensing for survey 
of earth resources and better weather -forecasting 
scientific space research etc. 

(iv) To provide for the changing satellite attitude so as to 

achieve its fixed orientation with respect to the moving 
line of sight appeared to be a major challenge involved 



40 


involvod in. tho use of the proposed nonsynchronous orbits 
for communications. The investigation showed that this 
challongu can he successfully net hy employing the passive 
attitude control method using the solar radiation pressure. 

5*2. Rcoommondations for Future Work 

There are several aspects related to the use of 24-hour 
non-synchronous satellites for cominunications which have not been 
considered in the present investigation. It would be worthwhile to 
carefully examine those aspects which may bo summarized as follows: 

(i) The ground station may require using more sophisticated 
system. It would, therefore, be essential to examine the 
effect of the resulting increase in the cost of ground 
component on overall economics. 

(ii) A transmission power control will become necessary to 
cuuntoract the variation of the range caused by orbital 
cllipticity. 

(iii) Several communication problem related to frequency coordina- 
tion, intormodulation are likely to be adversely affected. 
Besides, it would be desirable to extend the analysis to study 
the attitude control characteristics of the proposed solar 
model as affected hy the orbital inclination, roll and yaw 


motion 
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APffiNDIX 


Calculation of total velocity increment in transfer of satellites 
from parking orbit to 24-bour coplanar elliptic orbit. 

The geometry of satellite motion and its transfer from 
circular parking orbit to the final 24-hour elliptic orbit is shown 
in Fig. 2.1. It is easy to show that the circular velocity at P 
is given by 

T = 47^ (A.1) 

c 

The corresponding velocity at P along Hohmann transfer 
trajectory is given by 



The velocity increment required at p for injecting the 
satellite from parking orbit into the 24-hour elliptic orbit is 
therefore given by 


Ay^ = - 2a2/(a^+a^) - ^ a^/a^ 

(A. 3) 

On nott-diinensionalising the velocity increment, we get 
AVg = ’^2n2-1 - '^^2 




where 
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A .2 Calculation of total velocity increment in transfer of satellites 

from parking ortiit at an inclination with respect to the equatorial 
plane to f;oooynchronous orhit. 

As shovm earliest, velocity increment required at P to change 
the satellite orhit to Hohmann transfer trajectory can he written as 

AVp = '/\x /o-2_ Q/2a2/a^ - aa^/Ca^+a^) - ^ a^/a^H (A.5) 

Birailarly, tho velocity increment required at A in order to inject 
tho satollito into geosynchronous orhit can he obtained 

a^) >^+{2 sin i/2}^ 

(a. 6 ) 

where i is inclination of parking orbit, combine (A. 5) and (A. 6 ) 
results in the follov;ing non-dimensional expression for the total 
velocity increment 



(A. 7) 

where 

AVc =■ AVp / »'^/a2 
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Fig. 2-1 Geometry of 24-hour circular and elliptic orbits 














Fig. 3-5 System plots showing the effect of 
solar position angle on maximum 
librational amplitude 








55 



0 -Orbits 

Fig. 3 6 Typiccil r6spons0 showing tho 

affifude control through the proposed 

solar controller for (b) e = 0.1 
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Fig.3»7 System plots showing the attitude control 
characterstics as affected by the solar 
position angle. 
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Fig.3'8 System plots showing attitude control 
characterstics as affected by the 
eccentricity- 



= 0.05 
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Fig.3-9 System plots showing the effect of cnar.gir 
solar position angle from its nominal value 
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Fig.4-1 Librationai behaviour of satellite 
relative to the ground station as 
affected by eccentricity 




3tional behaviour of satellite relative to tl 
ons as affected bv solar narnmc?tf?rc: 




3 4-3 Librational behaviour of satellite 
relative to the ground as affected 
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Cj = C. i<P,Ki) 
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Fig. 4-4 Typical satellite response showing the 
control of librations of satellite about 
line of sight as achieved through the 
proposed solar controller 



Cj=Cj Ki) 
Cj = Cj (0 


Fig. 46 Librational control of satellites relative 
to the ground station as affected by 
position of the sun 



Fig.4-6 Librational control of satellite relative 
to the ground station as affected by 
the eccentricity 
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Fig.4-7 Typical response plots showing the relative 
attitude control performance 
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Fig. 4-8 System plots showing the effect of eccentricity and 
control gain ju on attitude control performance relative 
to the ground station 
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«IBFTC RKADE 

SUBROUTINE RK AU L ( X . XL . Y , Y 1 , 1 ) 


(, 

r 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 




« SUBROUTINE FQH SOLVING TFg SECOND ChCER « 

« ORDINARY DIFhfchENTlAL EGLATION-- 

* Y’ » = F( X, Y, Y' ) 

« X--AN ARRAY Uh X-AXIS VALLES « 

* XL — FINAL VALUE OF X « 

»- Y--AN ARRAY OF Y-AXIS VALLES 

« Y1--DY/DX « 

* Hl'-STEP LENGTH ■» 

» li— NO. OF lines already PRINTED AND « 

<» LEFT AS BLANK « 

« UW--INPUT-OUTHUT DEVICE FOR PRINTING « 

KK = 10 IF printing CF VaLLES IS NOT « 

* REQUIRED 

* nn-'-flag indicator « 

« =0 IF X AND Y in radians « 

« =1 IF X AND Y IN DEGREES « 

« =2 IF X AND Y IN RADIANS AND INTEGR- « 

« ATION IS FROM G TC 2«PI(OR NUlTIPLE « 

» OF PI ) « 

* FUNCTION F TO BE PRCVIDED BY USER « 

* ADAMS EXTRAPOLATION NETPCC USED TO SOLVE « 

* DIFFERENTIAL EQUATICN « 

* RUNGE KUTTA FOURTH CFDER METHOD USED FOR « 

» STARTING VALUES « 




REAL K1»K2,K3,K4 

DIMENSION FX(6) .DhLl(6).CEL2(6) ,BEL3^6) ,DELA(6) .D£L5(6) 
DIMENSION X(L) » Y(L) ,Y1CL) 

CCMMON/VKJ/Hl, i 1 »vjW,KH»NN 
PI =A .#ATAN(1 . ) 

PP = PI/180 . 

PI2=2.«PI 
H = H1 

IF (NN.EQ.l) 

H2 = Hi»H/2. 

H3=H/2. 

N =5 
M = N+1 
C 

C *tt« RUNGE-KUTTA METHOD ««« 

C 

IF(NN.EQ.l) Y=PP»Y 
DO 1 I=1*N 

IF(NN.EQ.l) X{i)=PP«X{I) 

FX( I)=F(X( I) ,Y( 1 ) »Y1( I ) ) 

K1=H2»F(X(I),Y(1)»Y1(I)) 

K2 = H2»F(XCI ) + H3,Y(U + H3«Yl(I)+Kl/4..Yin)+Ki/H) 

K3 = H2#F(X( I )+HS,Y( J )+h3«Y1 ( I ) +K1/ 4 , , Y1. ( I ) + K 2 /H ) 

K4sH2*F (X( I ) + H» Y( 1) +P«Y1( I )'+K3/ Y1 C 1 ) + K3/H3 ) 

IF(NN.EQ.l) X{l)=X(l)/FF 
K - I 1 

Y(K)sY( I)+H«Y1( i)+(Ka->K2'»K3)/3. 
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V 


Y1(K)=Y1( n+\K1+2.«K242.«H3 + K4)/3./H 

1 X(K)=X(1)+H1 
IF(NN,EQ.l) X (K )sPP»)( (K ) 

FX(K)sF(X(K), Y(K) ,Y1(K) ) 

IFCNN.EQ.l) X(K)=X<K)/FF 

C « «■» 

c ADAKS EXTRAPni. A! IOIS t'ETFCE 
C ««« 

2 DO 3 1=2, M 

3 UELKI ) = FX{I )-i-X(I-l ) 

DC 4 1=3, M 

4 DEL2( I ) = DEL1( I )-Dt:Ll( I-l) 

DC 5 1=4, M 

5 DEL3( I)=DEL2( I)-DEt2< I-n 
DO 6 1=5, M 

6 DE14( I)=DEL3( I)-DfcL3( 1-1) 

DEL5(M)=DEl.4(M)”DfcL4(K-l ) 

K=K + 1 

Y(K) = Y(K"l)+H»Yl(K-l)+F«F«(FX(H)/2.+nELlCI^)/6.+CEL?fM>/e.+ 
1DEL3( M)»19./18U . + DfcL4 (p ) « 3 , /32 , +DEL 5 ' '' ) «p 63 . / 10 OR 0 . ) 
Yl(K)=Yl(K-l)+H«(hXCF)+CELl(M)/2.+CbL2(M)*f;/12.^DEL3(P)«3./8.+ 
1DEL4(M)»251 ./72U.+CfcL5(P )«95, /2 88 . ) 

DO 7 1=1, N 

7 FX( I)sFX( I+l) 

X(K)=X(K-1) +Ha 
IFCNN.EQ.l) GO 10 8 
IF(NN.EQ.2) go !C 9 
FX(H)=F(X(K),YIK),Y1(K)) 

II-(X(K)-XL) 2,2,10 

a XX=X(K) 

X(K)=X(K)-FLOAT ah IX(X (K )/360, ) )«360 , 

IFCK'N.EQ.l) X(K)=PP«X(H) 

FX(H)=F(X(K),Y(K),Y1(X)) 

X(K )=XX 

U <X(K)-XL) 2,2,10 
9 XX=X(K) 

X(K)=X(K)-FL0A!(Ii-IX(X(K)/Pl2))«Pl2 

FX(P)=F(X(K).Y(K),Y1(K)) 

IF(X(K)-Xl.) 2,2,10 

10 IF(NN,NE.l) GO 10 12 
DO 11 1=1, L 

11 Y( I )=Y(I)/PP 

12 IF(KK.EQ,10) RhlUHN 
C «•»« 

C PRINTING the VALUES CF X AND Y 

C 

Jl=l 
J2 = 6 

13 11=11+1 

IF( II ,EQ.70) GO 10 16 
' WRITE(JW,15) {X(I),Y(I),I=w1,J2) 

Jl= Jl+6 
J2=J2+6 

IF(J2-L) 12,14,14 

14 02=02*6+1 
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WRITE(JW,15) (X(I),Y(n,ISw2,L) 

15 FCRYATC1X,6(F lU .1,H0 ,5 ) ) 

RETURN 

16 WRITE(JW,17) 

17 FORMATdHl,//) 

11=0 

GO TO 13 
ENT 


flIBFTC FAX 

SUBROUTINE MAX C Y , N , YF AX ) 

C fl SUBROUTINE USfcU TO FINE THE MAXIMUM « 

C « ABSOLUTE LIBRAIIONAL AMPLITUDE « 

C « Y--INPUT VECTOR « 

C » N--NO. OF ELEMtNTS U VECTOR Y « 

C « YMAX — MAXIMUM ABSCLITE CF VECTOR Y « 

DIMENSION Y(N) 

DO 1 1=1, N 
Z=A8S(Y( I ) ) 

1 IF(Z.GT.YMAX) YMAX=Z 
RETURN 
END 


«I8FTC FS 

FUNCTION F(THbiA,SI,Sll) 

C » FUNCTION SUB-PHOGRAM USED TO SOLVE THE « 

C * GIVEN DIFFEREMIAL EGLATICNS USING « 

C « CONTROL POLICIES FOR SATELLITE TO EXECUTE « 

C « ABOUT THE LOCAL VERTICAL « 

C * THtTA--POSITiON OF SATELLITE « 

c « s)--lI8Rational angle « 

C » S U--DERI VATI Vfc OF LIBRATICNAL ANGLE WITH « 

C » RESPECT TO ANGLE THETA « 

C * E--ECCENTRir. I i Y » 

C « KI--INERTIA PARAMETER 

C « PHI— -SUN POSIllON ANGLE « 

C » PI--CONSTANT 13.14159265) « 

C » AMUl, ANUl , AMU2. ANU2--CCNTRCLLER GAINS « 

C » CIMAX, C2MAX-~M AX IMUM PERMISSIBLE VALUE OF « 

C « SOLAR parameters « 

C » LL— FLAG INblUATOR TC LSE CNE CONTROL « 

c » STRTAGEY >1 A ONE TIME « 

C * C1,C2--S0LAR PAHAMETEFS « 

C 4nnnnte •»■»«■««««»**««* »*■»'****’** 

REAL KI 

COMMON E»KIiPh 1 ,Pi,AMll,ANUl,AMU2,ANU2,ClMAX»C2MAX»LL 

PI2=2,«PI 

A=l.+E 

8 = 1 dE«COS(THEl A) 

C= E«SIN(THETA) 
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E = THETA + SI-PHI + H1'<; 

D=D-FLOAT( IFIX(L/Pi2>)«Fl2 
IF(LL) 10,20,30 
10 C1=0, 

C2 = 0, 

GC TO 40 

20 Cl=-0,75<»PIfft»CQS(PHI)/A««3 
C2= 0 ,75«PI«fc«blN(PHI)/A««3 
GC TO 40 

30 Cl = -0 .75«PI*E<^CUS(PHI )« (E/A >«»3 
C2= 0,75«PI«E<fblN(PHI)«(E/A)»«3 
40 C11=-(AMU1*SI1+ ANU1«S I) 

IF(Cll-GT.ClMAX) C11=C1PAX 
IFCCll.LT. (-CIPAX) ) C11=-C1PAX 
IFCD.GT.PI .ANlJ.L.Lt.2.«FI) C11 = -C11 
C1=C1+C11 

C21= (AMU2«S11+ANU2«S1 ) 

IF(C21.GT,C2MAX) C21=C2PAX 
IF(C21.LT,(-C2MAX)) C21^-C2PAX 
lF(D.GT.PI/2, .ANU.D.LE.3,«PI/2. ) C21=-C2l 
C2=C2+C21 

G = (2.*C«(1.+SU )’3.«K I«SU(SI)<^COS(Sn)/B 

HsA«*3/8*»4«CCl«SlN(C)fl(AES(SlN(D) ) ) ”C2*CCS ^ D 5 «( APS( CCS ( E ) ) ) ) 

SI2=G+H 

F=SI2 

RETURN 

END 


«IBFrC FZ 

FUNCTION F(THfcl A,bI,SIl) 

C « FUNCTION SUB-PHOGRAP USEE TO SOLVE THE « 

C * DIFFERENTIAL hUUAlICNS USING DIFFERENT . « 

C « CONTROL POLICIhS FOR SATELLITE TO EXECUTE ® 

C « ABOUT THE LINE OF SIGHT « 

C « THETA--POSITION OF SATELLITE « 

c « si—librational angle « 

c * sii— derivative of lieraticnal angle with « 

c « RESPECT iO angle THETA « 

C « E--ECCENTRICi IT ^ 

C * KI--INERTIA PAKAMEfER « 

C » PHI--SUN POSIiiON ANGLE « 

C « PI--CONSTANT (3.14159265) « 

C » AMUl, ANUl, AMU2, ANU2--C'CNTRCLLER GAINS « 

C « C1MAX,C2MAX-”HAX INUH FERHISSI8LE VALLE OF « 

C « SOLAR FARAHETERS « 

C * AN--CONSTANI (6.628) « 

C « LL--FLAG indicator TC USE CNE CONTROL « 

C » STRTAGEY A! A ONE TIHE « 

C « C1,C2--S0LAR PARAMETERS « 

C ■» Wl-'-THET A-WE^f 1 * 

C « We--ANGULAR VELOCITY CF EARTH « 

c « T--TIME * 
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real KI 

COPHON E»KI ,HI . Af'Ll.ANLl, AMU2» AM;<,ClPA)i,CSPAX, AN.LL 
PI?=2,«PI 
A = 1 . +E 

8 = 1 .+E»COS( THE I A) 

C= E*SIN(THErA) 

D=THETA+SI-PHI 
D = n-FLOAT( IFlX(i;/PI2) )«FI2 
G = 1 . *E*E 

ai=an«g-i. 

IF(LL) 10 » 20, 30 
10 C1=Q, 

C2 = 0, 

GO TO 40 

20 Cl=(-0,75«PI«fc«COS(PH I ) + 2 . 25 *P I «K I / A1 »E*CCS C PH I ) w/s,-*3 
C2=( 0.75«PI«E»SlN{PI'n-2.25»PI«KI/AJ,-BEeSIN<PHI ) ^ ''A-n*3 
GO TO 40 
30 F=l,-E 

X= (SQRTC A)-SQR 1 (F) ) «S 1 N ( T F ET A ) / ( 2 . « ( SORT C A )«COS ( THE^ A/2 . )»»2+ 
lSQRT(F)»SIN(Thh lA/2. )a«2 ) ) 

Y=ARSIN(2.*X/(a .+X««2 ) ) 

W1=SQRT(G)/B*C+Y 

H-pH I -wl/Aj, 

Cl5(-.0.75»PI«EaCOS(H)*2.25«PI«KI/AlaEaCOS(l- ) )«(B/A>«a3 
C2= ( 0.75*PIaE*SIN(H)-2 .25«PI«KI/A1-»E«SIN (H ) )«(B/A)ae3 
40 C11=-(AMU1«SI1+ANU1«SI) 

IF(CII.GT.CIMaX) t;it=ClFAX 
IF(C11.LT.(-C1MAX)) C11=-C1PAX 
II- (C.GT.Pl.AND.L.Lfc.2.«FI) C11 = -C11 
C1=C1+C11 

C21= (AMU2«SI1 + ANU2«SI ) 

1F(C21.GT.C2MAX) C21=C2FAX 
IF (C21 .LT- (-C2MAX) ) C2l=-C2FAX 
IF(C.GT.PI/2. .AN!:.C.LE.3.«PI/2. ) C21=-C21 

C2=C2+C21 

G=(2.«C«(1.+SI1 )-3.«Kl«SlR(SI )*C0S(S1 > ) /B 

h = A<>»3/8»«4*(Cl*SIN (C)« { AES (SIN(D) ) ) -C2*C0S ( D ) « ( ABS* COS (D ) ) > ) 
S I2=G + H 
F=S I2 

return 

END 


«I8FTC 2IVAL 

SUBROUTINE Z I V AL ( T hET A 1 , Z I , N ) 

C «*«**««««« 


C a SUBROUTINE UShL 10 EVALUATE THE VALUE OF « 
C a ANGLE ZI FROM LiBRATICNAL ANGLE SI « 
C « THETAl--POSniUN Oh SATELLITE « 
C a ZI— INPUT VECIOR SI ANE RETURNS AS OUTPUT « 
C a VECTOR ZI « 
C a ALPHA--ANGLb bEIWfcEN THE LOCAL VERTICAL « 
C a and THh LINE OF SIGHT » 
C » N--NO. OF ELEMENTS IN VECTOR ZI a 





REAL KI 

DIMENSION THhTU(N),2ICN) 

COMMON E.KI ,PHC ,PI ,AMLl,AhLl,AMU2»ANU2»ClMAX,C2MAX. AN 

P12=2.«PI 

A=l.+E 

F=1 .-E 

G=1.-E*E 

DO 10 I=1,N 

THETAsTHETAK i ) 

B = 1 ,+E*C0S(THE I A ) 

C= fc*SIN(THETA) 

X = (SQRT(A)-SQR1 (F ))*SIN ( 7 HE TA ) / ( 2 . ■» (S CRT ( A ) «COS (THETA/2 .)««2+ 
lSGRT(F)«SlN(THt!A/2. )««2)) 

Y=ARSIN(2 .»X/ (1 . + X^i«2 )) +F 12 
W1=SQRT(G) /B«C+Y 
Wl=Wl-FLOAT(IFiXCWl/F 12 ))«F 12 
RN=AN*G/B 

ALPHA = ATAN(SIN( HI )/ (RN-CCS(^1 } ) ) 

10 2I( I)=2I( I)-ALHHA 
RETURN 
END 



